EPJ manuscript No. 

(will be inserted by the editor) 



Reweighting of the form factors in exclusive B —> X£is£ decays 

D. Cote, S. Brunet, P. Taras, B. Viaud 

Laboratoire Rene J. A. Levesque, Universite de Montreal, C. P. 6128, Succursale Centre-ville, Montreal, QC H3C 3J7, Canada 
May 27th, 2004 

Abstract. A form factor reweighting technique has been elaborated to permit relatively easy comparisons 
between different form factor models applied to exclusive B — > Xlvi decays. The software tool developped 
for this purpose is described. It can be used with any event generator, three of which were used in this 
work: ISGW2, PHSP and FLATQ2, a new powerful generator. The software tool allows an easy and reliable 
implementation of any form factor model. The tool has been fully validated with the ISGW2 form factor 
hypothesis. The results of our present studies indicate that the combined use of the FLATQ2 generator and 
the form factor reweighting tool should play a very important role in future exclusive \V u b\ measurements, 
with largely reduced errors. 

PACS. Form factors reweighting in exclusive semileptonic B decays 



1 Introduction 

Exclusive semileptonic B — > X u tvi decays can be used to 
measure the magnitude of the CKM matrix element V u b 
as their branching fractions (B.F.) are related to \V u b\ by 
the following relation: 

V X ■ TB 

where tb is the lifetime of the B meson and x is given 
by the relation : r th (B -> X u £ue) = x\V ub \ 2 , r th is the 
theoretical partial decay rate. 

The study of exclusive semileptonic B — * X u tvi decays 
offers some experimental advantages compared to an in- 
clusive study of all b — > ulvt decays, such as the possibility 
of keeping a higher fraction of the phase space and permit- 
ting a better background rejection. On the other hand, it 
deals with lower statistics and it is affected by large the- 
oretical uncertainties arising from the calculation of the 
form factors describing the strong interaction effects on 
the hadronization of the different X u final states. These 
uncertainties in the form factors lead to different predic- 
tions of the shape of the differential decay rate which, in 
turn, yield different predictions for the momentum spec- 
trum of the lepton £, the meson X u and its daughters. The 
subsequent varying efficiencies for the experimental cuts 
lead to uncertainties in the measured branching fractions. 
The theoretical uncertainties in the form factors also af- 
fect the values of a;. The resulting uncertainty in x is the 
largest source of uncertainty in the determination of \V u b\ 
from exclusive branching fraction measurements. 

In addition, in many analyses (e.g. those performed 
in BaBar), the simulated inclusive lepton spectrum of B 



decays does not agree with data. Since the B — ► Dlvf. and 
B — > D*£vi are the most abundant of all B decays, it is 
likely that at least part of the disagreement could arise 
from a wrong theoretical input to the simulation for these 
decays. Again, the most likely source of error comes from 
the form factors for these decays. Since the theoretical 
form factor predictions cover a rather large range, it is 
necessary to establish experimentally which theory best 
describes the data. This requires the possibility of varying 
the theoretical assumptions at the simulation level. 

A tool to be described in this paper has been created 
for this purpose. It will permit to switch easily between 
various form factor hypotheses and/or to vary the param- 
eters of a given hypothesis, within a full standard Geant4 
Monte Carlo (MC) simulation framework. The basic prin- 
ciple of this new tool is to generate events and run the full 
simulation and reconstruction sequence only once, with 
a given form factor hypothesis. Subsequently, the events 
thus generated are reweighted at the ntuple level with the 
values provided by a different form factor hypothesis. 

The relevant formulas and basic principles of the form 
factor reweighting technique are presented in Sect. In 
Sect. |3| the structure of the software tool developped for 
this purpose and how to use it are described. In that sec- 
tion, we also show how to incorporate new form factor 
hypotheses in the tool. Several histograms 1 that demon- 
strate that the tool is working properly are shown in Sect. 
01 A sample of kinematical distributions calculated with 
different form factor models are compared in Sect. [SJ The 
conclusions are given in Sect. El 

1 A far more extensive document is available from the 
authors. This document presents a large number of useful re- 
lations, part of the C ++ sofware tool developped in this work 
as well as many more histograms. 
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2 Technique for form factors reweighting 

2.1 Pseudo-scalar versus vector mesons 

The exclusive B — ► X u £ + vi decays 2 which could be stud- 
ied with this technique are: B + — ► ir° /rj/rj 1 / 'p° /u>£ + vi and 
B° — > 7r~ /p~£ + vi. As will be shown in Sect. EH the differ- 
ential decay rates are different for pseudo-scalar and vec- 
tor mesons. The B°, B^, ir^, it , r\ and rf mesons are all 
pseudo-scalar particles, while the p^ , p° and u> mesons are 
vector particles. As well, the most abundant B — > X c £ + vi 
decays, B — > D£ + vi and B — > D*£ + vg, involve D pseudo- 
scalar mesons and D* vector mesons. 




2.2 Kinematics of semileptonic decays 

A semileptonic B — > A^ decay is generally described 
by the following process. The B meson first decays into 
a virtual boson and an X meson which are emitted 
back to back in the B frame. The virtual W boson then 
decays to a lepton and a neutrino which are emitted back 
to back in the frame, while the X meson decays in 
various ways. The kinematics of such semileptonic decays 
can be completely described by three angles: 61, dy and \ 
defined in Fig. and by q 2 , the invariant mass squared of 
the virtual M /± boson. In terms of 4-momenta: 



Fig. 1. Definition of the angles 9i, 8v and % m the case of the 
B to vector semileptonic decay B — > D* tv, D* —* Dir. Ot is the 
helicity angle of the W ± boson, the angle between the direction 
of the W boson boosted in the B frame and the direction of 
the lepton boosted in the W ± frame (where the £ and v are 
emitted back to back). 9v is the helicity angle of the X meson, 
here the angle between the direction of the D* meson boosted 
in the B frame and the direction of the D meson boosted in 
the D* frame (where the D and ir mesons are emitted back to 
back), x is the angle between the W x I plane and the D* x D 
plane. 



q 2 = (pt + Vvf = (pb - Pxf 



(2) 



The four variables are totally uncorrelated. In the B frame, 
q 2 is also given by: 



q 2 = m 2 B + nix ~ ^msEx 



(3) 



where Ex is the total energy of the X meson. It is also of 
interest to note that the magnitude of the 3-momentum 
Px and q 2 are uniquely related by Eq. fj]in the B frame: 



Px = 



2.3 Form factors 



' (m 2 B + m 2 x — q 2 ) 2 



Ami 



•x 



(4) 



The matrix element of a semileptonic B 
can be written as UJ: 



X£f£ decay 



M(B -> X£v e 



-i-^V xb L»Hv 



(5) 



where Gf is the weak interaction's Fermi constant, V x b 
is either V u b or V c b depending on the final state meson, 
L M is the leptonic current and is the hadronic current. 
The leptonic current is well-known and can be calculated 
precisely using perturbation theory. 



2 Charge conjugation is implied throughout this paper, un- 
less explicitly stated otherwise. 



The hadronic current accounts for the strong interac- 
tions between quarks and gluons, and thus for the hadroniza- 
tion of the final state quarks into a A meson. The hadroniza- 
tion of a u quark and a d spectator quark into a -k~ meson 
in a B° — > tt^£ + V£ decay is a good example of such a sit- 
uation. In all these processes which involve the exchange 
of soft gluons, the strong interaction coupling constant, 
a s (p), is too large to allow the use of perturbative cal- 
culation techniques. Thus, even though the structure of 
hadronic currents is well known, such currents cannot be 
computed directly. However, they can be parametrized in 
terms of a small number of so-called universal Isgur-Wise 
functions, or form factors 0- To compute the form fac- 
tors, it is necessary to use either non-perturbative calcu- 
lation techniques such as lattice QCD [3], or approxima- 
tions of QCD. The approximation techniques can them- 
selves be split into two categories: 1) those, such as the 
Light Cone Sum Rules (LCSR) which are identical 

to QCD at some extreme limits but are a good approxi- 
mation of QCD in a restricted but known kinematic range 
and 2) those, such as ISGW2 which, instead of QCD, 
use approximate wave functions based on quark models 
for the mesons. The values of the form factors extracted 
from experimental data can then be confronted with var- 
ious models or used on their own. 

In terms of form factors, the hadronic current is given 
by a different expression depending on whether the B me- 
son decays to a pseudo-scalar or to a vector meson fi- 
nal state. In both cases however, the expression for the 
hadronic current can be simplified in the limit [S] of a 
massless lepton, assumed to be true when £ = e or p, (see 
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Sect. 14.41 for an estimate of the effect of this approxima- 
tion). 

For decays such as B° — > ir + £~Di where 7r + is a pseudo- 
scalar meson, the hadronic current is written, in the limit 
of a massless lepton as : 

ff" = (7r + (p')\urb\B0(p)) = f + (q 2 )(p + pY (6) 

where q — p — p' and f + (q 2 ) is the form factor describing 
the non-perturbative QCD effect. 

For decays such as B° — ► p + t~vi where p + is a vector 
meson, the hadronic current is written, in the limit of a 
massles lepton as 6 : 



H» = {p+(p',e)\u^(l- l5 )b\B°(p)) 



(7) 



m B + m p 
tub + rn p 



<P' a PpV{q 2 ) - (m B + m p )e^A 1 (q 2 ) 
(p + pYA 2 (q 2 ) 



an event. The total decay rate r does have a large effect 
PQ on the model dependency of B — ► Xlv analyses but is 
not used in the context of form factor reweighting. 

In the limit of a massless lepton, the differential decay 
rate of semileptonic B decays to a pseudo-scalar meson is 
given by [Hj: 



where Ai(q 2 ), A 2 (q 2 ) and V(q 2 ) are the three form factors 
describing the non-perturbative QCD effect. 

The equations for these two hadronic currents as well 
as their expressions without the massless lepton approx- 
imation are discussed in various papers, for example in 
Ref. 0. 

It is important to note that the leptonic current and 
the structure of the hadronic currents (Eqs. EJandQl fol- 
low directly from Lorentz invariance and are thus model- 
independent. The theoretical uncertainties in exclusive semi- 
leptonic decay analyses are only due to the uncertainties 
in the knowledge of the form factor(s): f + (q 2 ), Ai(q 2 ), 
A 2 (q 2 ) and V(q 2 ) which are model-dependent. These fac- 
tors nominally depend only on a single variable: q 2 . How- 
ever, in the case of vector meson decays, interference ef- 
fects between the Ai(q 2 ), A 2 (q 2 ) and V(q 2 ) form factors 
introduce an additional model dependence for the angular 
differential decay rates £TJ. 



2.4 Differential decay rates 

The differential and total decay rates are the observable 
manifestations of the underlying leptonic and hadronic 
currents. According to quantum field theory, the decay 
rate is given by a squared matrix element containing a 
combination of leptonic and hadronic currents. Just like 
for the structure of the current equations (Sect. the 
structure of the differential and total decay rate equations 
is also considered to be model- independent. Only the form 
factors appearing in the rates give rise to theoretical un- 
certainties in exclusive semi-leptonic decay analyses pTJ. 
These will be investigated. 

Our investigation will be greatly simplified by the use 
of a technique to reweight the form factors among the var- 
ious hypotheses under study. As will be shown in Sect. 12.51 
this reweighting is equivalent to a reweighting of the dif- 
ferential decay rates i.e. of the probabilities of generating 



dr(B -> Sl + vi) 
dq 2 d cos 9 1 d cos 9y d\ 



\V xb 



G 2 



I287T 4 



where q 2 



| Ps | 3 sin 2 9 e sm9 v \f+(q 2 
(8) 

and x have been defined in Sect. 12.21 ps 
is the 3-momentum of the final state pseudo-scalar meson 
in the B frame, f + (q 2 ) is the QCD form factor described 
in Sect. l2~3l and V x b is either V u t or V c b depending on the 
final state meson. It is often practical to use an expression 
where two or three of the angles are integrated out in 
which case Eq. |H1 becomes jjj: 



dr(B -> S£+ 



vi) 



dqd cos 9 1 
dr{B -► Sl+ 



= IK 



xb\ 



r 2 

32tt 3 



Ps| i W0 / |/- t V)r (9) 



V() 



dq 2 



= M, 



xb\ 



r 2 

24tt 3 



|psH/ + (<7 2 )| 2 



(10) 



Also, in the limit of a masless lepton, the differential 
decay rate of semileptonic B decays to a vector meson is 

® 



dr{B Vl + vi) 



114 



|2 3G||pvk 2 



dq 2 d cos id cos Oydx l ' J '"' 8(47r) 4 m 2 3 



(I - cos e ) 2 sin 2 v \H+(q 2 )\ 2 
+ ( I + cos 9 i ) 2 sin 2 6 v | H- (q 2 ) | 2 
+4 sin 2 9 1 cos 2 9 v \H (q 2 )\ 2 

— 4sin^(l — cos9e) sin#y cos^y cos xH+(q 2 )H (q 2 ) 
+4sin^(l + cos9g) sm9y cos9y cosxH_(q 2 )H (q 2 ) 
-2 sin 2 9 t sin 2 9 V cos 2 X H + (q 2 )H- (q 2 ) 

(11) 

where pv is now the 3-momentum of the final state vector 
meson in the B frame. The functions H + (q 2 ), H_(q 2 ) and 
Ho(q 2 ) are known as the helicity amplitudes of the vec- 
tor meson. They are related pQ to the QCD form factors 
A 1 (q 2 ), A 2 (q 2 ) and V(q 2 ) described in Sect. ESI Integrat- 
ing out the angles, Ea. rm becomes: 



x < 



dr(B -> Vl+vt) 



dq 2 



= Wx, 



> G 2 F |Pv|g 2 

QQrT^m 2 ^, 



{\H + ( q 2 )\ 2 + \H_(q 2 )\ 2 + \H (q 2 )\ 2 ) (12) 



As can be seen from Eqs.[S]and^] the form factors in- 
troduce a model dependence in the prediction of the shape 
of the differential decay rates of both pseudo-scalar and 
vector meson decays. This will be shown explicitly in Sect. 
With standard procedures, this theoretical uncertainty 
in the shape of the differential decay rate leads to a "the- 
oretical" uncertainty in the efficiency of the experimental 
cuts, and thus on the measured branching fraction. 
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2.5 Reweighting the probabilities of generating events 
among various form factor models 

A useful feature of the form factor reweighting tool is that 
events are generated, fully simulated and reconstructed 
only once, using a given form factor hypothesis. The prob- 
abilities of generating events are then reweighted to any 
other form factor model at the ntuple level. Such a tech- 
nique presents enormous advantages in terms of flexibility, 
time, CPU resources and disk space required, compared 
to generating, fully simulating and reconstructing sepa- 
rate data samples for each form factor model to be inves- 
tigated. 

Given the events generated with a certain probabil- 
ity by a specific form factor model G, the probabilities 
of generating events according to a different form factor 
hypothesis O, are obtained by applying a weight w to the 
probabilities of generating the events of type G. For a 
pseudo-scalar meson decay, the weights are defined as: 



dr(B^Slv) \ 
dq 2 d cos 9 e >° 



, dr(B->Slv) 
V dq 2 d cos 6 1 



(13) 



and for a vector meson decay as: 



w = 



, dr(B^Viu) s 
^ dq 2 d cos 6id cos 6yd\ ' 



dT(B^Vlv) 
dq 2 d cos 6 id cos 9yd\ ' 



(14) 



In this work, three different generators were used to 
calculate the initial probabilities for generating events. 
These probabilities were then reweighted to other form 
factor models. The generators are: 

— ISGW2: This generator, based on a quark model cal- 
culation 0, is extensively used in BaBar, Belle and 
CLEO. It is used for the simulation of generic BBbar 
events including that of the B — > X u tv decays. The 
differential decay rates are computed in this hypothe- 
sis using Eqs. 10 and 12 and the form factors given in 
Ref. (2 3 . 

— PHSP: This PHase SPace generator has been used for 
several decays in our work. It generates events with 
equal probability in all points of the phase space. In the 
context of this generator, the differential decay rate is 



dr(B^Xlv) _ 
dq 2 d cos 9t> d cos 9vd\ 



constant x 



given by the relation: 

|px|- This means that the generated cos 9i, cosOy and 
X distributions are flat while the differential decay rate 
decreases almost linearly with q 2 (see Eq. 0] and Fig. 

H. 

FLATQ2: This new generator pQ has recently been im- 
plemented in the BaBar software. It defines the prob- 
ability of each event as the probability given by the 
PHSP generator divided by the value of \px\ in the 
B frame, with a cut-off at |px| > 0.01 GeV/c. In this 
case, the differential decay rate is given by the rela- 



tion: 



dr{B^Xtu) 
dq 2 d cos Opd cos 9yd\ 



constant. This means that 



3 It should be mentionned that there are typographical errors 
in Ref. [2J. These have been corrected in our work. 



the generated cos 0i , cos 9v , X an d <? 2 distributions are 
all flat i.e. all the events are generated with an equal 
probablility throughout this 4-dimcnsional space. Of 
the three generators, this one is the most useful to ex- 
tract the form factors of B — > Xiv decays. In partic- 
ular, the fact that events are generated with an equal 
probability for the complete q 2 range is useful to eval- 
uate the efficiency of the experimental cuts, especially 
at high q 2 where most models predict very few events. 
The high q 2 events are of the utmost importance in 
the study of lattice QCD results. 

In cases where the initial distributions are generated 
with e.g. the ISWG2 generator, the distributions for any 
other form factor model, e.g. the LCSR one, will be ob- 
tained by applying, in the case of a B — > Stv decay, the 
following weights to the ISGW2 distributions: 



|T4 fc | 2 ^|p s | 3 sin 2 , 



\v xb 



\ftcSR(l 2 )\ 2 



|p s |3sin 2 ^|/+ GW/2 (g2 ) |2 



1/ 



ISGW2 



(15) 
(16) 



If the same B — » Siv decay distributions are generated 
initially with the FLATQ2 generator, the event-by-event 
weight has, in this case, a rather simple form: 



w 



sin 2 6e\p s \ 3 \ft CSR (q 



2\|2 



(17) 



3 The form factor reweighting software tool 
3.1 Outline 

The form factor reweighting software has been written in 
C ++ . The code is •practically self-contained, so that it can 
easily be used outside a specific framework (in that case, 
the CLHEP libraries need to be included). The code is 
written with an object oriented structure so that its dif- 
ferent sections are independent of each other, each section 
being a separate class. For example, the section for com- 
puting the kinematic variables of Sect. I2"2*l is independent 
of the one for computing the differential decay rate formu- 
las and weights which in turn is independent of the one for 
computing the form factors of the different models. This 
structure of separate classes having as input arguments 
either the LorentzVectors in the LAB frame or the output 
objects of the other classes in the reweighting software, 
controlled by simple user interfaces, yields a high degree 
of versatility. This design should allow an easy and reliable 
implementation of expected new form factor models. 



3.2 Software architecture 

The form factor reweighting software consists of two sepa- 
rate tools: XSLKin and XSLEvtFFWeight, built on three 
inheritance levels, each made up of one or more classes. 
The top level class of each tool is the user's interface from 
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which all the lower classes are inheriting, directly or in- 
directly. All variables or functions needed by the user are 
declared (known as "pure virtual functions" ) at this level 
and defined and used for computation at the lower lev- 
els. The software architecture is shown in Figs. |21 and |3J 
The three level structure of the XSLKin class diagram 
(Fig. [21 will probably be simplified in a future version of 
the sofware. On the other hand, the three level structure 
of the XSLEvtFFWeight tool (Fig. |3J) is very efficient in 
computing the reweighting among the various models to 
be investigated, as will be explained in Sect. 13.31 



3.2.1 The XSLKin tool 



The XSLKin tool computes the kinematic variables de- 
scribed in Sect. 12.21 namely: q 2 , 0£, 9y, x- The objects of 
this tool are used as input arguments by the XSLEvtF- 
FWeight tool to compute each event weight. The XSLKin 
tool can also be used as a standalone for other physics 
analyses. From the C ++ point of view, the tool's struc- 
ture is quite simple, as shown in Fig. [2] 



"interface" 
XSLKin 



XSLScalarKin 



q2(), thetaJO 



XSLVectorKin 



theta_v(), '■HQ-a . 



Fig. 2. XSLKin class diagram. The arrows mean an inheri- 
tance. 



"interface" 
XSLEvtFFWeight 



FromFLATQ2ToThisMo del() 

FromPHSPToThisModelO 

FromISGW2ToThisModelQ 




Fig. 3. XSLEvtFFWeight class diagram. The arrows mean 
an inheritance. 



clare the functions needed by the tool's user for reweight- 
ing an event. The second level of the diagram contains the 
classes that inherit from the XSLEvtFFWeight interface. 
In these classes, the functions declared in the interface are 
written out explicitly for the cases to be investigated i.e. 
for the B to pseudo-scalar meson decays in the XSLPseu- 
doScalarFF class and for the B to vector meson decays 
in the XSLVectorFF class. These functions are used only 
to calculate the weight required to reweight the proba- 
bility of generating an event from a given generator to a 
given form factor model. The form factors themselves as 
well as the kinematic parameters are computed in other 
classes. At the second level, the form factors are declared 
as pure virtual functions (GetFplusQ and GetAllFF()) to 
be denned and computed at the third level. The third 
level classes are then used to compute these GetFplus() 
or GetAllFF() virtual functions of the second level i.e. the 
f + (q 2 ) or Ai(q 2 ), A 2 (q 2 ) and V(q 2 ) form factor values as 
a function of q 2 for various models. 



The XSLKin user's interface is used to declare all the 
kinematic variables needed by users of this tool. The vari- 
ables are defined and computed at the second and third 
levels. The second level class, named XSLScalarKin, is 
used to compute the kinematic variables (q 2 and 8i) char- 
acterizing a B to pseudo-scalar meson decay. The third 
level class, named XSLVectorKin, is used to compute the 
two additional kinematic variables (Oy and x) needed to 
describe fully a B to vector meson decay. The XSLVec- 
torKin class already contains the values of q 2 and 9g in- 
herited from the XSLScalarKin class. 

This tool was created with an interface structure to 
allow the eventual addition of a second set of classes to 
compute the kinematic variables differently. It is not clear 
at present if this second set of classes will ever be used. 

3.2.2 The XSLEvtFFWeight tool 

The XSLEvtFFWeight tool architecture is shown in Fig. 
El The top level user's interface of this tool is used to de- 



3.3 How to implement new form factor models 

The XSLEvtFFWeight tool has been created with a three- 
level structure to make it efficient in computing the reweight- 
ing among various form factor models: the form factors 
are modified at the third level while the other two levels 
remain the same independently of the form factors used. 
The third level class inherits from either the XSLPseu- 
doScalarFF or XSLVectorFF class. The mathematical func- 
tions required to compute the form factor(s) as a function 
of q 2 are inserted in this class. 

4 Validation of the reweighting technique and 
its software tool 

In most simulations, the fully simulated events take into 
account the lepton's Final State Radiation (FSR). How- 
ever, the effect of FSR is not included in the reweight- 
ing formulas presented in Sect. l2.f)l As a consequence, the 
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Fig. 4. Events generated with the FLATQ2 (black), PHSP (pale grey) and ISGW2 (dark grey) generators: q 2 distributions 
for unweighted B — * r\lv decays (left panel) ; cos 8y distributions for unweighted B° —* plv decays (right panel) 



XSLEvtFFWeight tool could yield the wrong weights un- 
less corrected for FSR. We have found that the easiest 
way around this problem is to use the LorentzVector of 
the lepton calculated from the vectors of the other parti- 
cles of the decay as those are not affected by FSR. This 
solution turns out to be very effective even in cases where 
the FSR modifies the lepton's energy by several GeVs. As 
long as the same lepton is used to compute the kinematic 
angles and the event weight, the effect of the FSR on the 
form factor reweighting is negligible. 

The histograms shown in this section were calculated 
with such lepton Lorentz Vectors. They are extracted from 
1 million entries generated with no detector simulation. 
This generator information is what is used with full Monte 
Carlo simulation in real physics analyses. Only a small 
sample is presented in this section. More can be found in 
Ref. [I]. 

Our reweighting software has not been tested on any 
B — > X c lv decay but there is no reason to believe that 
the results will be any different from those obtained with 
the B — > X u tv decays. 

4.1 Properties of the generators FLATQ2, PHSP and 
ISGW2 

The properties of the three generators used in this work, 
and calculated with the XSLKin tool, are illustrated by 
the distributions displayed in Fig.Q] The distributions for 
all the decay modes of interest have been investigated ^ 
and found to be as expected. Fig. 21 (left panel) shows that, 
both, the PHSP and ISGW2 generators yield low statis- 
tics at high q 2 . This makes precise efficiency corrections 
difficult in this important region (particularly important 
for lattice QCD tests). Fig. 0] (right panel) shows that 



the PHSP and FLATQ2 cos Qy distributions are identical. 
This is also true for the cos^ and \ angle distributions 
PQ . Note that since the FLATQ2 generator yields flat dis- 
tributions for all the variables (9i, Oy, X, q 2 ) of interest, it 
allows a smooth reweighting over the full phase space. It is 
thus ideally suited to evaluate and correct the efficiencies 
as a function of these variables. To determine the analysis 
cuts required in realistic physics simulation, the FLATQ2 
and PHSP generated events have to be reweighted. 

4.2 Reweighting of B to pseudo-scalar meson decays 

Note that because of the three-level hierarchy of the XSLEvtF- 
FWeight tool, it is sufficient to fully validate the software 
with only one generator. For this purpose, the ISGW2 gen- 
erator was used. In this way, the XSLKin tool (Sect. l»T2~T|) 
and the classes XSLEvtFFWeight, XSLPseudoScalarFF, 
XSLV ectorF F, XSLScalarISGW2 and XSLVectorISGW2 
(Sect. 13.2.21) are all validated. In the implementation of 
any new form factor model, it is then sufficient to ensure 
the correctness of the new form factor equations at the 
third level of the XSLEvtFFWeight tool (see Sect. ETSj) . 

As can be seen in Fig.0 there is an excellent match be- 
tween the distributions generated directly with the ISGW2 
generator and those generated with the FLATQ2 or PHSP 
generators, and then reweighted to the ISGW2 form factor 
hypothesis with our software. 

4.3 Reweighting of B to vector meson decays 

Like for the reweighted pseudo-scalar meson decay re- 
sults, there is the same excellent match (Fig. EJ) when 
reweighting is applied to the distributions generated with 
the FLATQ2 or PHSP generators for vector meson decays. 
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4.4 Validity of the massless lepton approximation 



Comparing various distributions generated with a stan- 
dard ISGW2 generator to those generated with our form 
factor reweighting software allows us to test the validity of 
the massless approximation used in our reweighting soft- 
ware. This is the case since the standard code uses the ex- 
act formulas to compute the differential decay rates. Since 
the comparisons, done separately for electrons and muons, 
show the same good match, it can be concluded that the 
massless approximation is indeed justified for both, elec- 
trons and muons. An exemple of the quality of the match 
is displayed in Fig. 



1 16000^***^^ 



£ 14000* 

s 

■- 

~ 12000- 

c 

"° 

§ 10000 
8000 " 
6000^ 

4000: 

2000 1 



B n [iv 



ISGW2-EvtGen, no reweight 
FLATQ2 reweighted to ISGW2 



10 



15 



20 



25 



Fig. 7. Comparison of q distributions for B — > 717x1/ de- 
cays only. Events generated with a FLATQ2 generator and 
reweighted to the ISGW2 form factor hypothesis are shown in 
black while the unweighted events generated directly with a 
ISGW2 generator are shown in dark grey. 



5 Improvements in measurements of \V t 



ub 



With our reweighting software, it is now easy to investi- 
gate the predictions of various form factor models, and 
to evaluate their impact on the experimental study of 
B — > X u lv decays. So far, two form factor models for 
pseudo-scalar and vector meson decays have been fully 
implemented in our software: the ISGW2 model (used 
in the validation of the reweighting technique and its soft- 
ware tool) and a LCSR model from P. Ball et al. 0J[S]. 
Typical distributions for q 2 4 , cosOy and pe (Fig. EJ) de- 
duced from these two models display a large difference 
while the cos 9i distributions for pseudo-scalar meson de- 
cays (Fig. 03 left panel) are not model-dependent, but 
those for vector meson decays (Fig. |3l right panel) are. 
This significant model-dependence shows why the values 
of the branching fraction and of the CKM matrix element 
\V u b\, extracted from the study of exclusive semileptonic 
B meson decays, have such large theoretical errors. 



q is uniquely related to pxu 
Eq.El 



in the B frame as shown by 



The most precise value to date of |Vui,| obtained from 
the analysis of exclusive semileptonic B decays is the one 
given in Ref. 

\V ub \ = (3.27 ± 0.13 ± 0.19t°;^) x 1(T 3 

where the first error is statistical, the second systematic 
and the third theoretical. It is clear that the theoretical er- 
ror, of the order of 15%, due to form factor uncertainties, 
is by far the largest error. With our reweighting software 
tool, using our FLATQ2 generator, specially useful to de- 
termine the critical efficiencies of our experimental cuts 
as a function of q 2 , we aim to reduce the theoretical error 
to approximately 5%. The integrated luminosity required 
to achieve this goal depends on the technique used to re- 
construct the exclusive decay channel. With the neutrino 
reconstruction technique UJ, and the more than 250 fb~ x , 
collected at the T(4S) resonance in both BaBar and Belle, 
and of the order of 500 expected within the next two 
years, a smaller theoretical error should be possible in the 
near future. With the recoil techniques fTO], and for the 
other exclusive plv, ldIv, i]£v and r( Iv decay channels, a 
few additional hundreds of /6 _1 will be required to extract 
a similar theoretical precision on \ V u b\- 

With the soon to be available data from both BaBar 
and Belle, and with the neutrino reconstruction technique, 
a statistical error below 1% is possible. To achieve a the- 
oretical error at a few percent level will require tight con- 
straints on the form factors as a function of q 2 . This can 
be accomplished by two methods. 

In the first method, the measured differential decay 
distributions are compared to the corresponding ones pre- 
dicted by various form factor models, as shown e.g. in Fig. 
8. Our reweighting software tool will play a crucial role in 
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these comparisons since it will allow a quick and easy test 
of all models without having to regenerate a complete MC 
production for each model. In this method, the efficiencies 
of the analysis cuts are computed with the model used to 
generate the distributions. With very high statistics, the 
method should allow us to keep a single form factor model 
as being the one closest to reality, thereby reducing the 
theoretical error on \V u b\- 



In the second method, the form factors will be mea- 
sured directly with minimal model dependence. For exam- 
ple, for pseudo-scalar decays, the measured distributions 
will be fitted to Eq. 10 where the form factor is given e.g. 



EH by: 



n 



r-2 



1 



q 2 /m 2 (l 



j 2 /to 2 ) 2 



where r\ , r-2, and to are the parameters to be fitted. Other 
functions will have to be investigated to obtain the sys- 
tematic error engendered by the use of such functions. 
The small model dependence comes from the differential 
efficiencies, e(q 2 (,8e,8v,x))> of the analysis cuts used to 
generate the measured distributions. To evaluate these ef- 
ficiencies, the four kinematical variable distributions have 
to be binned, and the efficiencies are then determined for 
each bin. 

In the limit of infinitesimal bins, the efficiency e is inde- 
pendent of any model, as shown in Ref. [1]. But, of course, 
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in any practical analysis, with finite statistics and limited 
kinematical variable resolutions, the bin size will be finite. 
This introduces a small model dependency which must be 
taken into account. It can be limited by the appropriate 
choice of analysis cuts since it can be shown £Q that if 
the cuts are not correlated with the kinematic variables, 
then the efficiency is model-independent. The choice is 
greatly facilitated by the use of the FLATQ2 generator to 
determine the values of e, especially at very high values 
of q 2 . It is also facilitated by the use of the form factor 
reweighting tool to select the experimental cuts to limit 
the correlation and thus the model dependence. The gen- 
erator, combined with the software tool, will also yield a 
precise value of the small residual uncertainty due to the 
finite size of the bins. Once the parameters of the form 
factors are determined, they can be directly compared to 
any model predictions, thereby reducing the number of 
models and thus the theoretical uncertainty on \V u b\. 



9. BaBar collaboration, B. Aubert et al., Phys. Rev. Lett. 90, 
181801 (2003) 

10. D. del Re (private communication) 

11. P. Ball, R. Zwicky, |arXiv:hep-p h /0406232 vl, June 2004 



6 Conclusions 

A form factor reweighting technique and its software tool 
have been presented in this paper. Both have been vali- 
dated with the ISGW2 form factor model as illustrated 
by the excellent match between the distributions gener- 
ated with the FLATQ2 or PHSP generator, reweighted to 
the ISGW2 form factor hypothesis, and those generated 
directly with a ISGW2 generator. The object-oriented de- 
sign of the software tool allows an easy and reliable imple- 
mentation of any new form factor model, while optimizing 
the required CPU resources. The large differences, easily 
observed with our tool, in the distributions predicted by 
the ISGW2 and LCSR models for exclusive B -> X u lv 
decays show that a study of these decays will be very valu- 
able in extracting the values of the form factors. Our work 
leads us to expect that both, our novel FLATQ2 generator 
and form factor reweighting tool, will play a key role in 
the next generation of exclusive \ V u b\ measurements, with 
largely reduced errors. 
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